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We present new and exact "sign rules" for various spin-s anisotropic spin-lattice models. It is shown 
that, after a simple transformation which utilizes these sign rules, the ground-state wave function of the 
transformed Hamiltonian is positive-definite. Using these results exact statements for various expectation 
values of off-diagonal operators are presented, and transitions in the behavior of these expectation values 
are observed at particular values of the anisotropy. Furthermore, the effects of sign rules in variational 
calculations and quantum Monte Carlo calculations are considered. They are illustrated by a simple 
variational treatment of a one-dimensional anisotropic spin model. [PACS numbers: 75.10.Jm, 02.70.Lq, 
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Exact results for quantum many-body systems are 
rare. They are, therefore, always valuable as yardsticks 
against which numerical or semi-analytical methods may 
be measured. Such exact results acquire additional value 
when they can be used as actual input to improve one 
or more of the approximate methods. Such input is of- 
ten useful, for example, in tailoring a trial or starting 
many-body wave function to preserve certain exact prop- 
erties. In some instances they can even be practically 
essential for the successful implementation of particular 
techniques. An example is provided by the fixed-node 
quantum Monte Carlo (QMC) method, which provides 
a means to circumvent the infamous "minus sign prob- 
lem" inherent in simulating many-fermion systems, but 
at the cost of requiring good (and, ideally, exact) infor- 
mation on the nodal surface of the many-fermion wave 
function. In this paper we provide several exact results 
of such "sign rules" (which exactly define the nodal sur- 
faces) for some general classes of anisotropic spin-lattice 
models. These models are of considerable interest both 
in their own right and as models of real (low-dimensional) 
magnetic materials. 

In 1955 Marshall JjJ used a variational method to 
study the isotropic spin-half Hciscnberg antiferromagnet 
(HAF) specified by the Hamiltonian 
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where the sum on (i,j) counts each nearest-neighbor pair 
once and once only, and J is positive. Furthermore, the 
lattice was assumed to be bipartite (i.e., the lattice can be 
divided into two sublattices (A,B) such that all nearest 
neighbors of a site on one sublattice lie on the other and 
vice versa). The total number of atoms was TV and the 
spatial dimensionality was not restricted. In the course 
of this paper it was proven that the exact ground-state 
wave function can be written as 



I*) 



£<^> 



(2) 



where {\I)} are the usual basis states in the Ising repre- 
sentation. The coefficients {cj} were proven to have the 
property that, 
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where the coefficients {a/} are all positive real numbers 
or zero and </>(/) is the eigenvalue, with respect to the 
corresponding eigenstate \I), of the operator, 
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which counts all of the spin-half 'up' states on the A- 
sublattice. Lieb and Mattis M were able to prove that 
the ground state of the HAF was a singlet following the 
work of Marshall Q and Lieb et at. j§. 

Indeed, these results for the HAF model on a bipartite 
lattice are a consequence of a property of non-positive 
matrices. The Hamiltonian of Eq. (Q) is transformed in 
the following manner, 



H' = e^ffe"™* 
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It is seen that H' of Eq. (J5|) now contains only non- 
positive off-diagonal interactions with respect to the Ising 
basis states. Hence the ground state of the HAF corre- 
sponds to the eigenstate of largest magnitude eigenvalue 
which is positive-definite via the Perron-Frobenius the- 
orem M. Munro || and Parkinson fl were also able to 
extend these results to the spin-one biquadratic model 
on a bipartite lattice in various phases. They showed 
that the ground eigenstates are positive-definite and also 
non-degenerate. Finally, Klein M utilized this knowledge 
afforded by the Perron-Frobenius theorem for Heisenberg 
models to prove six theorems relating to ground-state fea- 
tures of these models. Amongst these theorems was one 
which stated that the ground-state expectation value of 
Si-Sj for these models is positive or negative depending on 
whether i and j are on the same or different sublattices, 
respectively. 

In this article, we present new sign rules for various 
spin-s anisotropic spin systems. The transformation of 



Eq. (|J) for these sign rules is utilized to force the off- 
diagonal interactions of the Hamiltonians to be always 
non-positive. A property of positive or negative semi- 
definite matrices is that the eigenvector corresponding 
to the eigenvalue of H' of largest magnitude (i.e., the 
ground-state energy here) contains only non-negative el- 
ements, and so the ground-state wave function is positive- 
definite. 

The first model that we consider is the XYZ model, 
which has a Hamiltonian given by 

hxyz = e{^ «t-j + Jit »r«? + ■% «?-?} . 

= El^ S i sZ 3 + bi 'j( S t S J + S i S t) 
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where 
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runs over all N sites 



h tim = Y, ^ lM 4A, +2E (4 + 4) ■ (8) 

[fci.fcz] i 



For the TIM, we place no restriction on the lattice type 
and k\ and k\ are allowed to run over all lattice sites with 
k\ =/= k,2- Furthermore, note that adding extra diagonal 
terms in Eqs. TO-H) (such as an external magnetic field 
or crystal field anisotropy in the z-direction), or letting 
i and j in Eqs. (o) and (H) run over the same sublattice 
but, in this case, explicitly restricting b iy j < in Eq. (H) 
and Jij < in Eq. ([?]), do not change the following sign 
rules. 

We now define the following operators, 
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hi = W& 

)/4. The index i in Eq. (g) 
on a bipartite lattice and j runs over all sites which are 
on the opposite sublattice to i. The square bracket [i,j] 
indicates that each bond is counted once and once only. 
The second such model that we shall consider is the 
single-ion anisotropy (SIA) model, given by 

#SIA = D 5>D 2 + E J ^ S « ■ S 3 > 

= E J « {**«? + l(sfsj + S r s +)| 

+ f E{^) 2+ ( 8 D 2 + 2(^+l)-(»f) 2 )}- (7) 

Again, the index i in Eq. (0) runs over all A^ sites on a 
bipartite lattice and j runs over all sites which are on the 
opposite sublattice to i. 

Finally, we also consider the transverse Ising model 
(TIM), which has a Hamiltonian given by 



where I a runs over all N/2 A-sublattice sites and / runs 
over all N sites on both sublattices. Note that for eigen- 
statc \I) = (g) ;=1 \si,mi), where sf\s u mi) = mi\si,n%i), 
m and itia have eigenvalues denoted by m(I) and 771^(7), 
respectively. Note that the following sign rules are fully 
defined by the form of the operator <fi which, in turn, 
gives the set of eigenvalues {</>(/)}. Furthermore, these 
eigenvalues determine the signs of the {c/} coefficients 
via Eq. (||) for the expansion of the ground-state wave 
functions, Eq. (0), of the anisotropic Hamiltonians of 
Eqs. (pl-13). The sign rules are thus defined by, 
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Furthermore, note that rule (B) is a reformulation, for 
general spin quantum number, of the Marshall-Peierls 
sign rule for the spin-half HAF given by Eqs. (|3|||). It 
is found that the first four sign rules apply to the XYZ 
and SIA models in the following regimes, 
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Rule 


XYZ model 


SIA model 


(A) 


hj < 0, dj < 


D < 0, J^ < 


(B) 


bij > 0, Cij > 


D < 0, Jij > 


(C) 


h,j < 0, dj > 


D > 0, Ji tj < 


(D) 


k,j > 0, a tj < 


D > 0, J-j > 
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for all values of i and j on opposite sublattices. (Note 
again that we do not allow the existence of any "frus- 
trated" interactions in Eqs. (|) or (|).) For the TIM, 
it is found that rule (A) applies when A < 0, and rule 
(E) applies when A > for all signs and strengths of the 
coefficients fik!,k 2 of the Ising interaction in Eq. dq). 

We once again note that the off-diagonal interactions 
of H' of Eq. (0) with respect to a complete set of Ising 
states are now always non-positive. Hence the ground- 
state eigenvector of H' contains only non-negative ele- 
ments and so may be written as |W) = ^2jdi\I), with 
a, > V I. 

It is also possible to prove exact relations (in some 
cases) for expectation values of various off-diagonal oper- 
ators. In order to illustrate this we define an XXZ model 
for a bipartite lattice in which we set J'f, — Jf 4 and 
J* a = 1 in Eq. (|6|). Hence we see from the table above 
that the sign rules are of type (A) for Jf : ■ < and of 
type (B) for Jf- > 0. We define the spin-spin correlation 
function in the ir-direction, given by 
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and the localized magnetization in the ^-direction to be, 



M x (l) = 
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For the XXZ model with Jf- > 0, it can be shown that if 
r in Eq. ( |l2| ) connects sites on opposite sublattices then 
G xx (r) < and if r connects sites on the same sublattice 
then G xx (r) > 0. Similarly for J X j > 0, it can be shown 
that the sign of M x (l) when I is a site on one particular 
sublattice must be the opposite to the sign of M x (l) when 
I is a site on the other sublattice. For J x j < 0, one can 
show that G xx {r) >0Vr and that M x (l) >0VI. For 
the transverse Ising model, it is found for any value of 
A that G xx (r) > V r. By contrast, it is found that 
M X {1) < V I for A > 0, and that M X {1) > V I for 
A < 0. We note that other such relations exist for the 
models presented here, although they are not discussed 
further in this article. 

In order to illustrate the usefulness of the sign rules 
given here, we now perform a simple variational calcula- 
tion for an anisotropic spin system which utilizes some 
of these sign rules. The system in question is a spin- 
half one-dimensional model J8| in which the anisotropy 
is in the x-direction such that Jf: = —A, Jf 4 — 1 , and 
Jf a = —1 in Eq. (|6j) and i and j are always nearest 
neighbors. Hence, from Eq. (Ill]) the sign rules are of 
type (A) for A > 1, of type (D) for -1 < A < 1, and 
of type (B) for A < — 1. The ground-state energy of this 
system using Eq. (13) may be written as, 
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In order to treat this model variationally we now em- 
ploy a Jastrow Ansatz, given by 



Cj = (-r/«<iin i + f(i,j)(PlP] 
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where the P^ and P^ are the usual projection operators 
of the spin-half 'up' and 'down' states respectively. The 
simplest approximation is to set the value of f(i,j) to 
be fi (a scalar variable) if i and j are nearest neighbors, 
and to be zero otherwise. The expectation value of Eq. 
(|l4J) may now be evaluated directly, although we explic- 
itly restrict the sums over all states in Eq. ( |l4| ) to those 
states in which m(I) is an even number in order to reflect 
the symmetries of the Hamiltonian ||. The variational 
ground-state energy is minimized with respect to f\ at 
each value of A, and the different sign rules for this sys- 
tem are utilized in separate calculations. We find that 
the best results for this variational Ansatz (for chains of 
length N < 16) are those which utilize the correct sign 
rule in each regime, and Fig. [j] illustrates these results 
compared to those of an exact diagonalization calcula- 
tion for a chain of length N = 12. (The accuracy of the 
variational calculations with N = 12 compared to ex- 
trapolation in the limit N — > oo via Pade approximants 
is estimated to be within about 4% for all A.) 
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FIG. 1. Results for the ground-state energy of an 
anisotropic spin-half quantum spin chain (described in text) 
of length N = 12. The variational calculation utilizes the 
Jastrow Ansatz of Eq. ( |15| ) in which the sign rules (A), (B), 
and (D) are used in separate calculations and these results are 
compared to those of an exact diagonalization of this system. 

This result is clearly true for any case in which one 
can prove an exact sign rule (e.g., those defined by Eqs. 
(ft) and (|l0|)) for a given model because any other choice 
would always contribute at least one non-zero positive 
contribution to the expectation energy from one of the 
'off-diagonal' contributions (i.e., I\ =/= I 2 ) in Eq. Jl4|). 
We may therefore infer that the best possible variational 
Ansatz for a given model must always utilize the correct 
sign rule in the correct regime. However, a direct evalu- 
ation of the sums in Eq. ( [14]) scales exponentially with 
increasing lattice size A, and so one must evaluate these 
summations using Monte Carlo (MC) techniques. In or- 
der to explain how this is achieved, we define a probability 
distribution for the Ising states {|/)}, given by 



P(I) 
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and the local energy of these states, given by 
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Equation (|l4J) may be equivalently written as, 
E = Y, P{I)E L {I) , 



(16) 
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where the sums over I\ in Eq. ( |l7| ) and / in Eq. (Hq) run 
over all, or possibly some symmetry-constrained subset 
(e.g., the subset of Ising basis states \I) for which m(I) is 
even), of the full set of Ising basis states. The MC proce- 
dure approximates the summation over / by performing 
a directed random walk with respect to \I) based on the 
probability distribution P(I). At each point of the MC 
'summation' procedure one 'steps' from one state \I) to 
another |/') (where |/') is one of a number of states ac- 
cessible to |/) via the off-diagonal elements of H) with a 



given probability which is dependent on P(I) and P{I')- 
Hence, it is possible to 'cover' all possible Ising states 
given enough run-time, although one samples the most 
important states the most often. Hence, in this manner 
an accurate approximation (to within statistical error) to 
the sum in Eq. ( Jig ) is built up. Also, the exponentially 
increasing problem with lattice size N is reduced to a 
problem that scales linearly with both TV and the num- 
ber of MC moves in a particular run. Hence, systems of 
larger lattice size may be treated. 

We note, however, that the sign rule does not affect the 
random walk which simulates the summation over \I) in 
Eq. (|18| ) in such a variational MC calculation because 
the probability distribution of Eq. (|16| ) is proportional 
to \cj\ 2 . By contrast, we may see from the expression 
for the local energy of Eq. ( |l7|) that each and every 'off- 
diagonal' contribution (i.e., / ^ I\) to the local energy 
can only be ensured always to be negative when the cor- 
rect sign rule in the correct regime is used. The average 
variational ground-state energy, evaluated using the MC 
procedure outlined above, is the average of the local ener- 
gies throughout the lifetime of the run. Hence the average 
variational MC ground-state energy can only ever be en- 
sured to be lowest (to within statistical error) when the 
correct sign rule is utilized within the variational Ansatz 
(e.g., in the Jastrow Ansatz of Eq. ([lq)). 

Let us finally turn our attention from variational MC 
estimates to full stochastic simulations of finite lattices by 
QMC techniques. All QMC methods at zero temperature 
basically project the exact ground state \ i $>) of a given 
many-body Hamiltonian H out of an initial trial state |$) 
not orthogonal to "J), by repeated applications of some 
suitable projection operator, G = G(H). This operation 
can be formally expressed as a path integral over many- 
particle trajectories in configuration space. Furthermore, 
it can thus be represented by a stochastic process, which 
may itself be simulated computationally by the random 
walks of a sampled set of independent "walkers" through 
the configuration space. The main limitation to the ap- 
plicability of such QMC techniques arises from the posi- 
tivity requirement, 



p(c,c/) = (c\G\c')^>0 



(19) 



on the probability p(c, c') for a walker to make a transi- 
tion from an initial configuration |c) to another configu- 
ration \c') along a random walk, where "5(c) is the wave 
function of the initial trial state. 

For fermionic systems, the minus sign problem asso- 
ciated with p(c, c') ensures that QMC calculations are 
much more difficult than for bosons. Thus, the computer 
time for bosonic simulations to give results of a specified 
accuracy scales algebraically with (i.e., as some positive 
power of) the system size, N. By contrast, when the 
wave function being sampled is not positive-definite, and 



hence cannot itself be regarded as a probability distribu- 
tion, no "exact" simulation method has been discovered 
that does not scale exponentially with N. Without some 
prior knowledge of the nodal surface a QMC simulation 
will always eventually relax to the corresponding bosonic 
wave function. For the typical spin-lattice problem for 
which no such exact knowledge is available most QMC 
attempts to alleviate the minus sign problem, such as 
transient estimation, are infeasible due to its severity. In 
these cases a recently proposed lattice variant M of the 
fixed-node QMC (FNQMC) method seems to pro- 
vide the only practicable QMC-based approach. 

Generally, results of FNQMC calculations are exact, 
within statistical errors limited solely by available com- 
puting power, only when the nodal surfaces of both the 
trial and exact ground-state wave functions coincide. 
Otherwise they yield variational upper bounds for the 
energy. On the other hand, other ground-state quantities 
such as the (magnetization or) staggered magnetization 
may be poorly estimated even when the energy is well 
estimated. How errors in the nodal surface affect such 
quantities is not well understood. In principle, however, 
their effects can be large. In this context the impor- 
tance of such sign rules as those presented here becomes 
clear, in that they permit a complete circumvention of the 
QMC minus sign problem. Indeed, it remains an open 
and challenging problem as to whether such complete or 
even partial sign rules can be found for other highly cor- 
related lattice systems which also suffer from the minus 
sign problem. Of particular interest, for example, would 
be to attempt to extend our work to both electron-lattice 
models and frustrated spin-lattice models. 
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